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Abstract We present results on the heating of the inter-cluster medium (ICM) by 
gravitational potential energy from in-falling satellites. We calculate the available 
excess energy of baryons once they are stripped from their satellite and added to the 
ICM of the hosting environment, this excess energy is a strong function of environ- 
ment and we find that it can exceed the contribution from AGNs or supernovae (SN) 
by up to two orders of magnitude in the densest environments/haloes. Cooling by 
radiative losses is in general fully compensated by gravitational heating in massive 
groups and clusters with hot gas temperature > 1 keV. The reason for the strong 
environment dependence is the continued infall of substructure onto dense environ- 
ments during their formation in contrast to field-like environments. We show that 
gravitational heating is able to reduce the number of too luminous galaxies in mod- 
els and to produce model luminosity functions in agreement with observations. 



1 Introduction 

Within the ACDM paradigm of structure formation, the growth of galaxies is gov- 
erned by dissipational processes such as cooling of gas and associated star forma- 
tion, as well as dissipationless processes such as mergers of stellar systems that are 
already in place e.g. l30l . The interplay of such processes regulates the properties of 
individual galaxies. At early times during the cosmic evolution (z > 2), when cool- 
ing times are short in low mass ifTUI and massive haloes J6] [18] dissipational pro- 
cesses associated with star formation either is disc-like systems or gas rich mergers 
ll20ll dominate the growth of galaxies. With the onset of downsizing in the star for- 
mation rate of galaxies at low redshifts e.g. l24l this picture changes, in particular 
for massive early-type galaxies that can accrete up to 80% of their stellar mass by 
z = from satellites fl4|. The importance of such 'dry' mergers for the growth of 
massive galaxies at late times has been investigated e.g. in ifTTl l22l [191 and recently 
linked to the size-evolution of massive early-type galaxies l27l [T5l l23ll . 

Feedback is generally assumed to play an important role in trends as the ones 
mentioned above e.g. l25l . In particular supernovae feedback has been considered 
in shaping the faint-end of the luminosity function |fT] [T6| . At the luminous end the 
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situation appears more complicated with several competing effects at play as e.g. 
AGN-feedback [3] or gravitational heating by in-falling sub-structure |[T7l l9l. 

The observed properties of the galaxy population in high-density environments 
show distinct properties in contrast to the field galaxy population. The population 
of spiral galaxies is gas-poor ETj and has smaller HI discs [28 1 compared to spiral 
galaxies in the field. Besides affecting the properties of individual galaxies, the mor- 
phological mix of the population as whole changes in high density environments as 
well. The so-called density-morphology relation shows and increasing trend in the 
fraction of early-type galaxies as a function local density [4j. The process of ram- 
pressure stripping (8) is able to remove gas from satellite galaxies orbiting in high 
density environments, and truncate disc sizes. The change in morphological fraction 
is generally attribute to a higher fraction of mergers during the formation epoch of 
the high density environment. 

While the above mentioned processes actively change individual galaxies we 
here address how the build-up of the environment passively affects the state of the 
central galaxy population living in it via the ability/inability to cool gas. Numerical 
simulations show that shock heating of gas falling into the cluster potential is able to 
heat gas initially to the virial temperature of the cluster 13. In addition conversion 
of gravitational potential energy of orbiting satellites is able to heat the inter-cluster 
medium (ICM) lTT7l l5l. In the following we will focus on the latter effect and include 
it into a semi-analytic model (SAM) to estimate its impact on the galaxy population. 



2 Model 

The results we present here are based on the SAM described in detail in |[T2l[T7l . 
using the three-year WMAP cosmology with Q = 0.27, Q A = 0.73 ,£2 b /£2 = 0.17, 
eg = 0.77, h = 0/71 [26). We briefly summarize here the main physical models that 
are important in high density environments, and that we have included in our SAM. 



2.1 Shock Heating 

Shock-heating of satellite gas takes place when a small satellite falls within the 
large potential well of a host dark matter halo and starts interacting with the hot 
halo gas. Generally SAMs assume that the shock heating is taking place instanta- 
neously, thereby removing all hot halo gas from satellites and adding it to the hot gas 
reservoir of the central galaxy 1101 . It has been that this process, if modeled instan- 
taneously might be too efficient resulting in too many passive satellite galaxies in 
high density environments. |29]. We adopt a simple prescription for the process of 
shock heating assuming that the rate at which mass is shock heated from the satellite 
is related to the halo dynamical time t^ yn viaM ~ M/f^ v „. This is effects allows the 
satellite to hold on to his hot gas for a longer period of time and to continue cooling. 
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Following [8), we assume that gas is stripped from the satellite once the dynamical 
pressure is able to over- come the gravitational force binding the gas to the satellite. 
In terms of energy deposited within the satellite gas, this can be approximated by 

Eram = Photon r\ (1) 

Here, we take \% ot as the velocity of the satellite perpendicular to its disc orientation 
and assume that the orbital velocity of the satellite is comparable to the sound speed 
Cgas of the hot gas. The efficiency of this process for a face-on disk should be max- 
imal, and the efficiency for an edge-on disk should be minimal; we take this into 
account by assuming that the disk orientation is random with respect to the infall 
direction. This is in good agreement with cosmological dark matter simulations that 
show that the spin vectors of merging dark matter halos are randomly aligned to 
each other fl"P3l . For simplicity, we calculate the density p/, or by taking the average 
density of hot gas within the host halos virial radius, and take r/, to be the character- 
istic half-mass radius of the gas within the satellite. The rate of stripped material is 
then calculated dividing E ram by the specific energy of the gas on which it is acting. 



2.3 Gravitational Heating 

The energy needed to strip baryonic material from satellites is provided by the con- 
version of gravitational potential energy, that the satellite gains during its infall on 
a slightly bound E^md ~ 0. Parabolic orbits have indeed been shown to be the most 
probable based on cosmological N-body simulations |[T3l . Gas that it stripped from a 
satellites can in principle have potential energy left, after it has been used up for the 
stripping process, and we assume here that this energy will be used to heat the inter- 
cluster medium. In practice we apply following equation to calculate the amount of 
gravitational potential energy left to heat the ICM: 

E g rav= £ hf gaS j(A<t>-AE stri p). (2) 

i'=l,n s <a 

The amount of stripped gas for each satellite ; orbiting within the host halo is M gU s,i 
and the total energy needed to strip the material from the satellites is AE str j p . The 
heating of the ICM based on Eq.|2] predicts that for the same sum of stripped ma- 
terial from massive satellites, heating is less efficient than if that material would be 
stripped from less massive satellites. 
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Fig. 1 The integrated efficiency of gravitational heating e ? ,„ r as a function of host halo mass. 
Galaxies living the the most massive haloes always experience the largest contribution from grav- 
itational heating, due to the continued infall from satellites. Shaded areas show the expected range 
in efficiencies for SN and AGN-feedback. 



3 Results 

The time integral of Eq.[2]over the evolution of a high density environment can be 
quite substantial and, if expressed in terms of E grav ^ ot — e grllv m t c 2 , we find values 
for 

Egrav ranging from a few time 10 8 to a few times 10 4 in galaxies with stellar 
masses ~ 10 10 and ~5 x 10 12 M Q , respectively. To emphasize that this trend is 
driven by the environment, we show in Fig. [T]the trend with halo mass which we 
take as a proxy for the environment. Massive haloes have more gravitational heating 
during their evolution than less massive ones. One can understand this behavior by 
considering the accretion of satellites onto halos. In general, the most massive halos 
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Fig. 2 Top: the probability distribution for the ratio of heating to cooling rates as function of 
environment. The x-axis shows the temperature of the hot halo gas as a proxy for environment. 
Bottom: The average gravitational heating rate as a function of environment. 



at any redshift had the largest accretion rates in the past, which explains the large 
amount of gravitational heating and its dependence on halo mass. It is worth compar- 
ing gravitational heating to other common heating mechanisms such as supernovae 
and AGNs. Comparing e grm , to e SN ~ 2.8 x 1(T 6 and eBH ~ 6.5 x -6.5 x 1(T 7 
shows that in general, gravitational heating is more efficient than supernova feed- 
back only in galaxies larger than a few times 10 11 M© and in halos more massive 
than 5 x 10 12 M ol j ot at z = 0. This regime corresponds to massive field galaxies and 
extends into group-like environments. For even more massive galaxies and dark halo 
masses larger than 10 13 M Q , gravitational heating starts dominating over proposed 
AGN feedback rates. 
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To illustrate the contribution from gravitational heating, we display the contours 
of the conditional probability for the ratio of heating to cooling that individual galax- 
ies experience in a given environment. The top panel in Fig.|2]shows the probability 
contours for central galaxies at z = 0. 1 in our simulations. We translate the deposited 
energy per unit time into a heating rate, labeled heat in Fig. reffig2, by calculating 
the amount of cold gas that can be heated to the virial temperature of the dark matter 
halo in which the galaxy resides. The cooling rate, labelled cool is the standard ra- 
diative cooling rate that we calculate in our SAM based on the prescription in fPHl . 
Figure [2] shows clearly that gravitational heating dominates in environments with 
temperatures above ^keV. The heating rate for the central galaxies can be up to 10" 
times larger than the cooling rate, and in the most dense environments the heating 
rate becomes 10 4 M yr _1 . The heating rate shows a clear environmental depen- 
dence that reflects the higher abundance of satellites that contribute to gravitational 
heating. From these results, one expects that star formation in central galaxies of 
dense environments will be shut down. 

To further illustrate the effect of gravitational heating on the galaxy population 
we show in Fig. [3] the luminosity function of galaxies at z ~ 0. We find a good 
match to the observed one by 0. The effect of gravitational heating is strongest for 
massive galaxies. These tend to live predominantly in high density environments, 
subject to gravitational heating. For comparison we also show a model without the 
contribution of gravitational heating. This model severely overproduces the number 
density of luminous galaxies. 



4 Conclusions 

We here show the effects that the surplus of gravitational potential energy from in- 
falling baryonic matter can have on the heating of the ICM. Galaxies in low density 
environments are generally not affected by gravitational heating as the contribution 
is small compared to other feedback sources like SN or AGN feedback. We find 
that in the most dense environments, like massive groups or clusters, the constant 
infall of substructure is able to provide enough energy to level, and even exceed the 
losses due to radiative cooling, hence stopping any cooling in such environments. 
The integrated contribution from gravitational heating is in cluster environments 
clearly larger than the contribution from AGNs suggesting that this could be an im- 
portant feedback source. This conclusion is supported by looking at the luminosity 
function of galaxies. Including gravitational heating is able to reduce the amount of 
too luminous galaxies and to bring them in agreement with observations. To further 
test the importance of gravitational heating high-resolution simulations of galaxy 
clusters are required. Such simulations should be able to reveal self-consistently the 
importance of the effects discussed here. 
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Fig. 3 The r— band luminosity function of galaxies in our model including gravitational heating 
(filled circles) and without (dotted line), compared to the SDSS results (solid line) by |21 
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